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ABSTRACT 

The Planck satellite, successfully launched on May 14th 2009 to measure with unprecedented accu- 
racy the primary Cosmic Microwave Background (CMB) anisotropics, is operating as expected. The 
Standard Model of the Universe ( "concordance" model) provides the current realistic context to ana- 
lyze the CMB and other cosmological/astrophysical data, inflation in the early Universe being part of 
it. The Planck performance for the crucial primordial parameter r, the tensor-to-scalar ratio related 
to primordial B mode polarization, will depend on the quality of data analysis and interpretation. 
The Ginzburg-Landau approach to inflation allows to take high benefit of the CMB data. The fourth 
degree double well inflaton potential gives an excellent fit to the current CMB+LSS data. We evaluate 
the Planck precision to the recovery of cosmological parameters, taking into account a reasonable toy 
model for residuals of systematic effects of instrumental and astrophysical origin based on publicly 
available information. We use and test two relevant models: the ACDMr model, i.e. the standard 
ACDM model augmented by r, and the ACDMrT model, where the scalar spectral index, n s , and r are 
related through the theoretical "banana-shaped" curve r = r(n s ) coming from the Ginzburg-Landau 
theory with double-well inflaton potential. In the latter case, the analytical expressions for n s and r 
are imposed as a hard constraint in the Monte Carlo Markov Chain (MCMC) data analysis. 

We consider two Cg— likelihoods (with and without B modes) and take into account the white noise 
sensitivity of Planck (LFI and HFI) in the 70, 100 and 143 GHz channels as well as the residuals from 
systematics errors and foregrounds. We also consider a cumulative channel of the three ones. We 
produce the sky (mock data) for the CMB multipoles Cf T , Cj E \ Cf E and Cf B from the ACDMr 
and ACDMrT models and obtain the cosmological parameter marginalized likelihood distributions 
for the two models. Foreground residuals turn to affect only the cosmological parameters sensitive 
to the B modes. As expected, the likelihood r distribution is much clearly peaked near the fiducial 
value (r = 0.0427) in the ACDMrT model than in the ACDMr model. The best value for r in the 
presence of residuals turns to be about r ~ 0.04 for both the ACDMr and the ACDMrT models. The 
ACDMrT model turns to be very stable, its distributions do not change by including residuals and 
the B modes. For r we find 0.028 < r < 0.116 at 95 % CL with the best value r = 0.04. We also 
compute the B mode detection probability by the most sensitive HFI- 143 channel. At the level of 
foreground residual equal to 30% of our toy model only a 68% CL (one sigma) detection is very likely. 
For a 95% CL detection (two sigmas) the level of foreground residual should be reduced to 10% or 
lower of the adopted toy model. The lower bounds (and most probable value) we infer for r support 
the searching of CMB B mode polarization in the current data as well as the planned CMB missions 
oriented to B polarization. 

Subject headings: cosmology: cosmic microwave background - cosmological parameters - inflation; 
methods: data analysis; space vehicles. 
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1. INTRODUCTION AND WORK OUTLINE 

The Planck satellite^ was successfully launched on May 14th 2009 to measure the primary Cosmic Microwave Back- 
ground (CMB) temperature and polarization aniso tropics on the whole sky with unprecedented accuracy. It is now in 
normal operation, with the expected performa nces dPlanck Collaboration 2006: Bersanelli et aJ.J[2010t lMandoles i~et al.l 
120101: lLamarre et all 120101: iMaffei et all I2010D . Planck will improve the measurement of most cosmological parame- 
ters by several factors with respect to current experiments, in particular the Wilkinson Microwave Anisotropy Probe 
(WMAP) satellited. The expected CMB polarization measurements from Planck will allow to push both (E and B) 
polarization results well beyond the present knowledge and considerably constrain the tensor (B modes) to scalar 
ratio parameter r, if not to obtain a detection on it. In this respect, the way of extracting and physically interpreting 
cosmological parameters (once the CMB data cleaned from the different astrophysical foregrounds) will be important. 
In other words, the Planck actual performance for the crucial primordial parameter r will depend on the adopted 
physical modeling and on the quality of data analysis and interpretation. It is then important and timely to make 
forecasts for the Planck determination of r and other cosmological parameters taking into account the theoretical 
progress in the field and WMAP results. 

The Standard Model of the Universe (or "concordance" model) provides the current realistic context to analyze 
the CMB and other cosmological/astrophysical data. Inflation (quasi-exponential accelerated expansion) of the early 
Universe is a part of this model and one important goal of CMB experiments is probing the physics of it. Inflation solves 
the shortcomings of the decelerated expanding cosmology (horizon problem, flatness, entropy of the Universe), and 
explains the observed CMB anisotropics providing the mechanism for the generation of scalar and tensor perturbations 
seeding the large scale structures (LSS) and primordial (still undetected) gravitational waves {B mode polarization). 

The current CMB + LSS data support the standard inflationary predictions of a nearly spatially flat Universe with 
adiabatic and nearly scal e invariant initial de nsity perturbations. These data are validating the single field slow- 
roll inflationary scenario (iKomatsu et al.l 120091) . Single field slow-roll mo dels provide an appealing, simple and fairly 
generic description of inflation ( Dodelsonl 120031 iBovanovskv et ahl I2009T ) . The inflationary scenario is implemented 
using a scalar field, the inflaton with a potential V(tp), sclf-consistcntly co upled to the space-time metric. In the 
effective theory based on the Ginzburg-Landau (G-L) approach to inflation l)Bovanovskv et aT1l2009D . the potential 
is a polynomial in the field starting by a constant term. Linear terms can always be eliminated by a constant 
shift of the inflaton field. The mass (quadratic) term can have a positive or a negative sign associated to unbroken 
symmetry (chaotic inflation) or to broken symmetry (new inflation), r espectively. The fourth degree double-well 
inflaton potential gives an excellent fit o f the present CMB + LSS data (jBovanovskv et aT1l2009D . A cubic term does 
not improve the fit and can be omitted (iDestri et al.ll2008aD . Adding higher order terms with additional parameters 
does not improve significantly the fits (jDestri et al.ll2009|) . The G-L framework is not just a class of physically well 
motivated inflaton potentials, among them the double and single well potentials. This approach provides the effective 
theory for inflation, with powerful gain in the physical insight and analysis of the data. The present set of data with 
the effective theory of inflation favor the double well potential (jBovanovskv et aLll2009t IDestri et al.ll2008af ). Analyzing 
the present data without the relation between r and n s does not allow to discriminate among different classes of models 
for the inflaton potential in the considered framework. Although the G-L effective theory approach to inflation is quite 
general, it predicts prec ise order of magnitude estimates for n s , r and the running of the spectral index dn s /dlnk 
(jBovanovskv et al.ll2009h 



here TV ~ 60 is the number of efolds since the cosmologically relevant modes exit the horizon till inflation ends. The 
WMAP values for n s and the upper bounds for r and dn s /dhik agree with these estimates. Since in this framework 
the estimated running, dn s /dlnk ~3x 10~ 4 , is very small, in this paper we will concentrate on n s and r. 

In this work, we evaluate the accuracy in the recovery of the cosmological parameters expected from the Planck data. 
First, we do this forecast without including the systematic effects of instrumental and/or astrophysical origin or their 
coupling, affecting the Planck measurements, and then by including the systematic effects. In this study we exploit 
the Planck sensitivity and resolution at its three favorite cosmological channels, i.e. at the frequencies of 70, 100, and 
143 GHz. Table Q] reports the Planck perf ormance at these frequencies, based onlPlanck Collaboration] (120061) and, for 
the LFI channel at 70 GHz, as updated in lMandolesi et al.l (|2010MBersanelli et all (j2010P ). lSandri et al.l (|2010h . These 
sensitivities do not include the degradation in accuracy that could come from various sources of systematic effects, of 
both instrumental and/or astrophysical origin, or their coupling. In Sect. [5] we discuss the current published estimates 
for the residuals of systematic effects and foregrounds affecting the Planck CMB measurements: straylight, main beam 
asymmetry, leakage, time constants, glitches, and foregrounds. In general, we do not use in this work a precise (still 
not completely available) description of the considered systematic effects, but only suitable representations of them, 
as described in Sect. [H This is done in a parametric approach, identifying the corresponding levels at which the 
control of the systematic effects is necessary not to spoil the Planck data scientific accuracy. We technically implement 
this rescaling with a multiplicative constant on the residuals of the systematic effects on the CMB multipoles Ct. 
Obviously, the real analysis of Planck data will have to properly consider all possible systematic effects of optical, 

6 http://www.rssd.esa.int/planck 

7 http://lambda.gsfc.nasa.gov/ 
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thermal, and instrumental (radiometric and bolomctric) origin, with an even better accuracy than those achieved in 
past projects. In parallel, a significantly improved separation of CMB from astrophysical components will be needed, 
a task in principle possible for Planck thanks to its wide frequency coverage. 

Instrumen t al sys tematics on CMB tensors-to-scalar have been studied by iHu et all (|2003|) ; iShimon et al.l (|2008|) ; 
lYadav et all IpOlTO . 

We use and test two relevant models: the ACDMr model, that is the standard ACDM model augmented by the 
tensor-to-scalar ratio r, and the ACDMrT model, that is the ACDMr model in which the double-well inflaton 
potential (see Eq. ((TJ) in the next section) is imposed. Namely, n s and r are constrained by the analytic relation 
r = r(n s ) to lay on the theoretical banana-shaped curve (the upper border of the banana-shaped region Fig. [1}. 
The novelty in the MCMC analysis of the CMB data with the ACDMrT model is in the fact that we impos e the 
analytical expressions for n s and r derived from the inflaton potential as a hard constraint (|Destri et a l. 2008aj). We 
take both models, ACDMr and ACDMrT, as fiducial models in our Monte Carlo Markov Chains (MCMC) simulations 
to produce the corresponding skies (mock data) . In the ACDMr model the independent cosmological parameters are 
flbh 2 , Vl c h 2 , 9, t, A s , n s and r, while all other independent parameters are assumed to vanish, e.g. Q, v = 0, or have 
the standard values, e.g. w = — 1. The aforementioned ACDMrT model includes the same parameters but with n s 
and r not being independent, but related by the curve r = r(n s ) as widely discussed in Sect. El We produce one sky 
(mock data) for the anisotropy CMB multipoles Cj T , Cf E , Cf E and C BB from the ACDMr model and from the 
ACDMrT model, with the parameters in Table [2] We describe the detailed procedure in Sect. [5] We run Monte Carlo 
Markov Chains from this sky and obtain the marginalized likelihood distributions for the cosmological parameters 
(fifth 2 , £l c h 2 , 9 r, Oa; Age of the Universe, z re , H , A s , n s and r) in the two test models ACDMr and ACDMrT . 
We study the independent ACDMr parameters with the mock data produced from ACDM (first row of Table [2]) and 
the independent parameters of both ACDMr and ACDMrT with the mock data produced from ACDMrT (second row 
in Tabled). The fiducial values, r = 0.0427 and n s = 0.9614 correspond to the best fit to the CMB-LSS data with the 
ACDMrT model using the double-well inflaton potential expressed by Eq. (QJ. 

Namely, these are the best fit values to r and n s within the Ginsburg-Landau effective theory approach. Not using 
the Ginsburg- Landau approach, lower bounds for r are n ot obtained and the best fit value for r can be much smaller 
than r = 0.04 (jKinnev et al.ll2008t [Peiris fc Eastherll2008l ). 

We consider two choices for the Cg— likelihood, one without the B modes and one with the B modes and take into 
account the white noise sensitivity of Planck (LFI and HFI) in the 70, 100 and 143 GHz channels (|Planck Collaboration! 
120061 ). We also consider a cumulative channel whose x 2 is the sum of the x 2 ' s of the three channels above. When 
using different channels in the MCMC analysis, we use different noise realizations while keeping the same sky, that is 
the same realization of the Gaussian process that generated the primordial fluctuations. In our MCMC analysis we 
always take standard flat priors for the cosmological parameters. In particular we assume the flat priors < r < 0.2 in 
the ACDMr model and < r < 8/60, where 8/60 ~ 0.133 is the theoretical upper limit for r in the ACDMrT model. 

We performed the MCMC simulations using the publicly available CosmoMC codc0 (jLewis & Bridle 2002) interfaced 
to the Boltzmann code CAME0 (see lLewis et al.l <|2000f ) and references therein). 

Our findings without including the systematic effects are summarized in Figs. EMS where the marginalized likelihood 
distributions of the cosmological parameters are plotted for several different setups. In Tables [3] to [5] we list the 
corresponding relevant numerical values. Clearly, in the case of the ratio r, due to the specific form of its likelihood 
distribution, it is more interesting to exhibit upper and lower bounds rather than mean values and standard deviations 
as in Tables [3] and [4] We report the upper bounds and, when present, the lower bounds in Tables [5] and [6] Our 
conclusions without including the systematic effects are: 

• The upper bound on r and the best value of n s do not require to include the B modes in the likelihood, and can 
be obtained with the ACDMr model alone, (i.e. r < 0.068 and n s = 0.9549 at 95 % CL ). See Tables El [5] and 
Fig. [3] The inclusion of B modes, for a non vanishing fiducial value, (r = 0.0427), allows peaked marginalized 
distributions for r and a lower bound for r. Sec Table [6] and Figs. El 111 and [5] We obtain 0.013 < r < 0.045 at 
95 % CL in the ACDMr model, with the best values r = 0.0240, n s = 0.9597. This shows a substantial progress 
in the f orecasted bounds for r with respect to the WMAP+LSS data set for which r < 0.20 in the pure ACDMr 
model (jKomatsu et al.|[200l [2010h . 

• Lower bounds on r and most probable r values arc always obtained (with or without the B modes) with the 
ACDMrT model. See Tables El H M and Fig. [SJ The ACDMrT model provide well peaked distributions for r on 
nonzero values r ~ 0.04. We obtain r > 0.039 at 68 % CL and r > 0.030 at 95 % CL in the ACDMrT model. 

In Sect. [7] we include in the forecasts the systematic effects discussed in Sects. El and 0] Our conclusions including 
the systematic effects and foreground residuals are: 

• The likelihood distributions with and without B modes result almost the same when including the residuals. 
Only the cosmological parameters sensitive to the B modes appear to be affected by the residuals, namely, r, z re 
and r. The main numbers are displayed in Tables [5] and El 

8 http:/ /cosmologist.info/cosmomc/ 

9 http://camb.info/ 
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• The marginalized likelihood r distribution for fiducial ratio r = 0.0427 is much clearly peaked on a value of r 
near the fiducial one in the ACDMrT model than in the ACDMr model (compare Figs. [7] and [8]). In any case, 
the best value for r in the presence of residuals is about r ~ 0.04 (near the fiducial value) both for the ACDMr 
and the ACDMrT models. The ACDMrT model turns to be robust, it is very stable (its distributions do not 
change) with respect to the inclusion of residuals (and they do not change neither with respect to the inclusion 
of B modes). The main numbers are included in Tables [5] and El With the ACDMrT model we have for r at 95 
% CL: 

0.028 < r < 0.116 with the best values r = 0.04 n s = 0.9608 . 

It must be stressed that, in the ACDMrT model, future improvements in the precision S on the measured value of 
n s alone will immediately give an improvement dr/dn s S on the prediction for r as well as for its lower bound. Better 
measurements for n s will thus improve the prediction on r from the T, TE and E modes even if a secure detection of 
B modes will be still lacking. 

In order to assess the probability for Planck to detect r we also compute the B mode detection probability by the 
most sensitive HFI-143 channel; this is done in Sect. 17.21 We extract 10 5 skies obtaining the corresponding multipolcs 
A[ m from the ACDMrT model according to the procedure described in Sect. [5j adopting r = 0.0427 as fiducial value. 
We compute all the corresponding likelihood profiles only for r and their interesting properties, like the most likely 
value r maxi the mean value r mean , the standard deviation Ar max of the r max distributions, the skewness and the 
kurtosis, (which measures the departure from a Gaussian likelihood), Fig. |H1 We finally compute the 99% CL, 95% 
CL, and 68% CL lower bounds for r. The probabilities of detection of r are displayed in Fig. [TO] At the level of 
foreground residual equal to 30% of the considered toy model, only a 68% CL (one sigma) detection is very likely. 
For a 95% CL detection (two sigmas) the level of foreground residual should be reduced to 10% of the considered toy 
model, or lower. 

Lensing acts on the B-modes as a contamination by transforming E-modes into B-modcs. It is a frequency indepen- 
dent effect while residuals are frequency dependent. Lensing weakens the signal around I ~ 90 where the primordial 
B-modes peak but not in the small £ modes range where the reionization bump dominates. On the other hand fore- 
ground residuals are larger at small I than at I ~ 90. Namely, residuals and lensing affect the detection of B-modes in 
complementary ways, with the effect of residuals stronger than that of lensing. As a consequence, lensing plus residuals 
can spoil the detection of r even when residuals are assumed at the 30% level of the considered toy model. On the 
contrary, lensing in the absence of residuals still allows a detection of r. For example, several MCMC simulations 
show that our lower bounds on r are not significantly affected by lensing in the absence of foreground residuals. Let 
us make clear, at any rate, that lensing was not considered in the analysis of the r— detection probability in Sect. 17.21 
Finally, it should be clear that if the theoretical constraint r = r(n s ) of the ACDMrT model is imposed on the MCMC 
analysis, r has always well defined lower bounds regardless of lensing and/or residuals. 

The forecasted probability of detecting r is based on the statistics of the shape of the r -likelihood. This shape 
determines whether a detection of r can be claimed with a given confidence level. But real CMB experiments can 
observe only one sample: the observed sky. So, the possibility of inferring r from one single (albeit very large) sample 
depends on the sample itself, and therefore, whether r will be or will be not detected depends also of a question on 
luck. 

In addition, the results for many skies presented in Sect. !7.2l show the consistency of our whole approach to determine 

r. 

Finally, in Sect. 17.31 we consider the bias effect in the foreground residuals implemented as a linear perturbation 
affecting the C[s and explore how the cosmological parameter distributions are affected by the bias. We implement 
two extreme cases: in case (i) the bias fluctuates randomly around zero and in case (ii) the bias fluctuates around a 
non-zero value, staying significantly non-zero. In case (i) the cosmological parameters are practically unaffected while 
in case (ii) the peaks of the cosmological parameter distributions are shifted within one or two sigmas of the WMAP 
values. In particular, r is not anymore detected in case (ii). 

The best and mean values reported here for r and the other cosmological parameters do not correspond to the true 
sky data but to mock skies generated from the MCMC simulations as explained above. Nevertheless, the deviations 
between the best and the fiducial values are relevant indicators for r as well as the lower and upper bounds and the 
standard deviation. The fact that the fiducial and mean values of r are very close and that Ar max coincides with the 
mean value of the standard deviation of r indicate that Planck can provide detections of high quality. 

More in general, our results support the quest for B mode polarizat ion in the current CMB data and future B 
oriented polarization missions under study by both ESA0 and NAS^EB (|de Bernardis et all [20091 : iBock et al.ll2006f) . 

2. FITTING CURRENT CMB + LSS DATA WITH THE GINZBURG-LANDAU EFFECTIVE THEORY OF INFLATION 

As discussed in the introduction, the effective theory of inflation within the G-L approach gives precise order of 
magnitude estimates for the spectral index n. s , t he ratio of tensor to scalar fluctuations r and the running of the 
spectral index dn s /dlnk {Boyanovsky ct al. 200§). 



10 http://www.b-pol.org/index.php 

11 http://cmbpol.uchicago.edu/ 
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Within the contex t of the G-L effective t heory of inflation, the wor k in lBovanovskv et al.l (2006) , D estri et al.l (|2008af) , 
IDestri et all (|2008bD . IDestri et all (|2009T >. iBovanovskv et all (|2009f ) showed that: 

• The small inflaton selfcoupling arises naturally as the ratio of the inflation energy scale and the Planck energy. 
The inflaton mass is small compared with the inflation energy scale. 

• The amplitude of the CMB anisotropics sets the energy scale of inflation to be M ~ 10 16 GcV for all generic 
slow-roll inflationary potentials. 

• Double-well inflaton potentials give the best fit to CMB+LSS data. Basically, the inflaton potential must have 
a negative second derivative at horizon exit which favours double-well potentials over single well potentials. 

• For double-well quartic inflaton potentials, the best value for the ratio tensor to scalar fluctuations is r ~ 0.05 
with the lower bound r > 0.023 (95%) CL in the case of the quartic double-well potential. The novelty in the 
MCMC analysis of the CMB+LSS data that leads to these results is in the fact t hat we imposed the analytical 
expressions for n s and r derived from the inflaton potential as a hard constraint (|Destri et al.ll2008al) . 

• Higher order double- well inflaton potentials arc investigated in IDestri et al.l (|2009( ). All r = r(n s ) curves for 
double-well even potentials of high order fall inside a universal "banana-shaped" region B, Fig. [TJ 

The fourth order binomial potential provides the simplest double-well potential best reproducing the CMB+LSS 
data within the G-L effective theory approach: 

/ x A / 2 m 2 \ 2 , y ( M \ 4 M 2 

y ( „ ) = _(^__) , A=^(_), m __ (1) 

where tp is the inflaton field, A stands for the quartic coupling, [y being the corresponding coupling of order one), m 
is the inflaton mass. Notice that the quartic coupling A is proportional to the ratio M/Mpi to power four and hence 

very small as stated above for all the inflaton self-couplings. 

Adding higher order terms with additional parameters does not improve significantly the fits (|Destri et al.l 120091 ) . 

In IDestri et al.l (|2009h it is found that the r = r(n s ) curves for double-well inflaton potentials in the G-L spirit fall 
inside the universal banana region B depicted in Fig. [T] 

The lower border of the universal region B is particularly relevant since it gives a lower bound for r for each 
observationally allowed value of n s . For example, the best value n s = 0.964 implies from Fig. [T]that r > 0.02L The 
upper border of the universal region B tells us the upper bound r < 0.053 for n s = 0.964. Therefore, we have within 
the large class of potentials inside the region B 

0.021 < r < 0.053 for n s = 0.964 . 

Moreover, the fourth order double-well potential represented by Eq. ([!]) is the simplest and G-L stable inflaton 
potential reproducing very well the present CMB+LSS data. 

Not using the Ginsb urg-Landau approach, the lower bounds for r arc not obtained. iKinnev et al.l (|2008ft ; 
iPeiris fc Eastherl (|2008[) do not use the Ginsburg-Landau approach, do not find lower bounds for r and cannot 
exclude values for r much smaller than r = 0.0427. 

It must be noticed that our present analysis shows that values r <C 0.0427, (and hence very small B modes) are 
outside the possibilities of detection by Planck. 

Future improvements on the precise value of n s alone will immediately give an improvement on the theoretical 
prediction for r as well as for its lower bound. An improvement S on the precision of n s implies an improvement 

dr 
dn s 

on the precision of r. According to IDestri et al.l ()2009f ) . at n s = 0.964 from r = r(n s ) we have: 
dv 

— — = 4.9 on the upper border of B (fourth degree double well) and 
dn s 

dv 

— — = 1.35 on the lower border of B. 

dn. 



Better values for n s will thus improve the prediction on r from the T, TE and E modes while a secure detection of B 
modes is still lacking. 



6 



Burigana et al. 



3. PLANCK SENSITIVITY 



In this study we exploit the Planck sensitivity and resolution at its three favorite cosmological channels, i.e. at the 
frequencies of 70, 100, and 143 GHz. 
Table [1] reports the Planck performance at these frequencies, based on iPlanck Collaboration! ( 20061) but consistent 



with the most rec ent pre-launch measurements of the HFI channels at 100 GHz and 143 G Hz ( Lamarre et al J 120101 : 
Maffei et al1l2010D and, for the LFI channel at 70 GHz, as updated in iMandolesi et all ()2010D . lBersanelli et all (|2010l) . 



Sandri et alJ (j2010f ). Notice that the LFI sensitivity reported here includes also the fluctuations of the 4K reference 



load, since it is obtained on ground-based calibration performed under realistic conditions. The resolution at the 
various frequencies comes from accurate optical simulations. Also notice that these numbers arc likely conservative, 
i.e. in principle a further refinement of tuning could return into an improvement of in-flight sensitivity. Almost four 
surveys have been adopted in this work. 

4. RESIDUALS FROM SYSTEMATIC EFFECTS AND FOREGROUNDS: TOY MODEL 

The sensitivities presented above do not include the degradation in accuracy that could come from various sources 
of systematic effects, of both instrumental and/or astrophysical origin. 

In this section, we discuss the current estimates publicly available for the systematic effects affecting Planck mea- 
surements. 

4.1. Straylight 

Planck achieves very g ood side lobe rejection thanks to its telescope design (jSandri et al.ll200"H 120101 : IMaffei et al.1 
120101 : iTauber et alJl2010T ). In spite of this, the main source of contamination at large angular scales or at low multipolcs 
comes from the so-called straylight effect, i.e. the signal entering in the lobes at various angular distances from the 
main beam. It can be distinguished in straylight from the intermediate beam, i.e. at angular distance of few degrees 
from the main beam, and from the far beam, i.e. at angular distance from the main beam larger than some degrees. 
The main sources of straylight are the Galactic emission and the CMB dipole. 

The straylight from the intermediate beam introduces a sort of smearing of signal around that observed by the main 
beam. Detailed studies show that it is important only close to the Galactic plane, a region typically excluded from 
scientific analysis through suitable masks, while it is significantly less of the straylight from the far sidelobes in all the 
other sky regions, and, ultimately, for the recovery of the CMB angular power spectrum. We will then consider only 
far sidelobes in the following estimates. 

Notice that, if the optical behaviour is well known it is possible to subtract to high precision this effect from the 
data by simply evaluating it on the observed sky by means of convolution codes taking into account the effective 
observational strategy, with the only (small) limitation introduced by the receiver noise. In practice, this correction 
is limited by the accuracy in the knowledge of optical behaviour. We will assume for numerical estimates an effective 
uncertainty of ~ 30% in the beam response in the sidelobes, implying that the amplitude of the spurious effect 
remaining in the data is about one order of magnitude smaller than the original effect. 

Therefore, assuming the (conservative) sidelobc levels computed for Planck at frequencies of 70 and 100 GHz, a 
reasonable esti mate for the residual strayl i ght from the Gala xy, rescaled from the computations carried out for the 
original effect (|Burigana et al.ll2001al feOoVlSandri et al.ll2010f ). is 

C t - ^ + ^ - 8 x 10~ 4 iiK 2 for I < 10 , (2) 

2 7T 

and 

C t iil±ll „ 2.5 x 10~ 4 fiK 2 for I > 11 . (3) 

2 7T 

Notice that at frequencies v > 70 GHz dust emission is the main Galactic diffuse foreground, while at 30 and 44 
GHz also free-free and synchrotron emission are relevant. Of course, the straylight effect is larger at 30 and 44 GHz 
but we could neglect them in this study for cosmological parameter estimation. Assuming similar sidelobe levels at 
143 GHz, given the typical dust frequency behaviour at millimeter wavelengths T„ almost proportional to v 2 , we have: 

C t i£±ll „ 3.2 x 10- 3 [iK 2 (4) 

2 7T 

Ct ~ 1 x 10- 3 »K 2 (5) 

A 7T 



for £ < 10, and 



for £ > 11. 

Th e other relevant source of straylight contamination is the CMB dipole ([Burigana et al.l 120061 : iGruppuso et al.l 
|2007|) . Note that, for symmetry reasons, this effect is significant only at even multipoles while it is negligible in 
practice at odd multipoles. Again, rcscaling the results obtained for the original effect, we derive a suitable range for 
the estimate of the residual contamination 

£ (£ + 1) 

d - 0.016 ~ 0.16 [iK 2 (6) 

2 7T 
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for even multipoles (with a typical value of 0.048 fJ.K 2 ) and 

«^ = ,7) 

Z 7T 

for odd multipoles. The larger values apply at lower multipoles (up to about 10), the lower ones at higher multipoles. 
Although the exact value depends on the particular considered receiver, we assume these estimates constant with 
frequency, as in the case of sidelobe levels approximately constant with frequency, being the CMB signal frequency 
independent (in equivalent thermodynamic temperature). 

Notice that, at the frequencies of 70-^143 GHz, for even multipoles dipole straylight is larger than Galactic straylight. 

In general, the contamination from straylight can be modeled to first approximation as an additional spurious excess 
of power. In principle, one could also include perturbations multipole by multipolc (or multipolc band by multipolc 
band) of the above estimates to avoid a modeling in terms of a simple analytical form for the spurious additional 
power. 

The straylight effect in polarization mainly depends on the (non-perf ect) balance of the str aylight in total intensity 
in the coupled receivers used to extract the Q and U Stokes parameters (jBurigana et a l. 2004). On the basis of optical 
simulations, we could assume relative differences of few tens per cent, which should reduce of about one order of 
magnitude the original effect with respect to that in total intensity. On the other hand, the modeling and verification 
of optics in polarization is much more complex than in total intensity. We then expect that it will be more difficult 
to use optical predictions to subtract this effect into the data and assume a residual effect similar to the original one. 
Therefore, we estimate that the amplitude of the residual effect in polarization will be similar to that in total intensity. 

4.2. Main beam asymmetry, leakage, time constants, and glitches 

Another potential systematic effect that needs to be kept under control is the effect of the antenna beam profile in 
the main beam. Beam profiles exhibit a deviation from perfect circular symmetry, in the ran ge from a few percent up 
to ~ 30% in the case of the lowest frequency channels (|Sandri et al.ll20i"ot iMaffei et al.ll2010t) . Main beam distortions 
are in principle a source of concern as they can bias the estimated power spectra in the high I regime, and hence 
affect the likelihood models and cosmological parameters. This happens for two reasons. PlancKs scanning strategy 
is not isotropic but has a preferred direction, roughly coincident with ecliptic meridians. In the first place, this fact 
makes the i space equivalent window function of the beam rather non trivial, an d difficult to estimate analytically 
(though approximate analytical solutions do exist, see e.g. iFosalba et al.l (|2002j) ). Secondly, and more important, 
Planck estimates the Stokes linear polarization parameters by combining measurements taken from different detectors. 
The beam asymmetry renders the contribution to the intensity I unbalanced even when the same pi xel is observed, 
becaus e of the different orientation and shape of the beams. In turn, this can create I to Q, U leakage (jAshdown et al.l 
l2009a| ) and produce biases in the polarization power spectra. 

Fortunately, the beam pr ofiles for Planck h&ve be e n measured very well d u ring ground testing c ampa igns and will be 
cross checked in flight (see iBurigana et all (|2001bD . iNaselskv et al.l (j2007f ). iHuffenberger et all (|2010f ) and references 
therein). Furthermore, analytical and semi-analytical machinery exist to estimate the I equivalent window function 
asymmetric beams, once a beam profile is known and a scanning strategy assumed. These methods compute, for each 
multipole I the beam coupling matrix betwe en all power spectra (t hus taking leakage into account), starting from an 
approximate model of the scanning strateg y (lAshdown et afll2009af) that can be refined performing signal-only Monte 
Carlo simulations for the CMB component (jAshdown et aljr2009b[ ). While a thorough analysis of the accuracy of these 
procedures have not been performed yet, it is fair to expec t that main beam distortions will not be a major source of 
systematic contamination for Planck (see e.g. iRocha et al.l (j2009h ). 

Not even a satellite experiment as Planck can safely assume to use the entire sky for CMB analysis. Incomplete sky 
coverage can induce leakage of the E polarization modes into B modes if a sub-optimal power spectrum estimator is 
employed. While this effect is not connected to the beam, but rather of geometrical origin, it is worth mentioning here 
because it may trigger spurious detecti on of B modes. In a realistic analysis, the leakage effect is corrected from the 
beginning by using pseudo Ci methods (|Hivon et al.ll2002l ) that are the standard choice for power spectrum estimation 
in the high I (> 30). Pseudo Ce methods correct for leakage by means of coupling kernels: in parti cular, the E 
and B mode pseudo-spectra exhibit correlations that need to be accounted for (see, e.g., appendix A in iKogut et al.l 
(|2003h ). At low multipoles, pixel based methods are normally used to compute directly the likelihood function with- 
out assuming power spectrum estim ation as an intermediate step. Pixel based methods do not suffer from leakage 
(jTegmark fc de Oliveira-Costal[200l . 

The bolometric detectors of HFI exhibit a non trivial transfer function, that distorts the signal both in amplitude and 
in phase. Qualitatively, the effect on amplitude is akin to a first order lo w pass filter arising from the detector's intrinsic 
time constant modified by electro-thermal effects (|Lamarre et al.ll2010D . Knowledge of the filter function allows one to 
deconvolve out its effect on timelines at the price of increasing (slightly) and distorting the high frequency noise leve l, 
which becomes non white ('colored noise'). These measurements will be performed in flight (jLamarre et al.l [20101 ). 
Any residual error would have an impact similar to beam smearing along the scan direction, so it contributes to beam 
asymmetry which, as stated above, can be accounted for with high confidence. 

Another potential source of concern for bolometric detectors is that they are sensible to cosm ic ray hits, that creat e 
glitches in the timeline, i.e. they are always seen as positive spikes in the bolometer signal (jLamarre eFaD[20i3), 
followed by a tail also due to the bolometer's time constant. These events can be detected and flagged on the timelines. 
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Residual effects (due to undetected glitches) can be kept under control (|Masi et al.ll2010D for the cosmological analysis 
by relying on angula r power spectra obtained cross correlating different detectors (also known as cross-spectra see e.g. 
iPolenta et all pOOl ). 

4.3. Residuals from foregrounds 

The most important source of contamination in a CMB experiment like Planck will come from the residual of 
astrophysical foregrounds. In fact, although the wide Planck frequency coverage, possibly complemented by WMAP 
maps and ground-based and balloon-borne experiment data, is particularly advantageous for a precise removal of 
astrophysical signals from the maps and the accurate mapping of CMB anisotropics, nevertheless we expect that a 
certain level of residual contamination will remain into CMB maps, particularly in polarization. Many methods of 
component separation, each with its own pros and co ns, have been and are continuously elaborated for the analysis of 
Planck multifrequency maps (see iLeach et al.1 ()2008f ) and references therein). In the present work, we are interested 
in the residuals from astrophysical foregrounds affecting the recovery of the CMB angular power spectrum. It is 
typically given as a difference between the input CMB angular power spectrum and the CMB angular power spectrum 
estimated after the component separation layer. In general, it is not so meaningful to provide a description of foreground 
residuals different at different frequencies, since, by definition, the component separation layer exploits exactly the 
multifrequency mapping of the sky. Thus, the estimate adopted in this work have to apply to the whole set of frequency 
channels. 

Different methods show different residuals at various ranges of multipoles. The multipolc dependence, or, in other 
words, the shape of this residual also depends on the consider ed method. 

Concerning residuals for the T mode, recent simulations (|Leach et al.l 120081) show residual shapes only slightly 
dependent on the multipolc, with amplitudes in the range 

C,«±»~¥±>* I X>, (8) 
Zir Zir 

the exact value depending on the method and, for each method, on the particular multiple band, with typical variations 
of about 30-40%. In this work, we model this T spurious power as flat in C( £ (£ + 1). 

Galactic polarized foreground (mainly from diffuse synchrotron and dust emission) affects CMB angular power spec- 
trum recovery more significantly in polarization than in temperature. We expect that their residual after component 
separation will take partially memory of the original shape of the foreground power spectrum, in particular at large 
scales where they show much more power than the CMB. Again, different methods give different residuals, regarding 
both multipolc dependence and amplitude. In this work, we model the foreground residual for the TE mode and E 
and B modes (assumed to be equal, B — E) as the sum of two shapes, the first one (dominant at low multipoles) 
described by the foreground shape properly rescaled in amplitude, the second one constructed from the foreground 
shape properly rescaled in amplitude and changed in slope. 

Fig. [2] displays our "starting" conservative models for the residuals in TE mode and in polarization modes. We 
let us the freedom to simply rescale them with multiplicative factors in order to address typical level of foreground 
residuals for which the impact on our cosmological aim is not critical (see also Sect. I4.5j) . 



4.4. Additional noise versus bias 

All the systematic effects discussed above, coming from instrumental effect, sky signal, or from their coupling, can 
be considered in two different schemes. 

In the first, simplest case, they can be treated as sources of spurious additional noise power, i.e. they do not 
introduce a bias affecting the recovery of the estimation of the CMB angular power spectrum but they increase our 
uncertainty in its recovery. Therefore, the effect can be modeled adding in quadrature the quoted Cg of the power of 
the residual systematics to those coming from sensitivity, resolution, and cosmic plus sampling variance. This approach 
is equivalent to assume that we will be able to properly model and subtract a correct estimation for the systematic 
effects so that only a statistical uncertainty in their subtraction will affect the data. 

In another, more critical approach, one can assume to miss the correct estimation of the spurious effects. Their 
systematic effects will be then much more dramatic, i.e. they will introduce also a bias in the estimation of the CMB 
angular power spectrum. This case can be modeled "perturbing" the Ce to be compared with the exact model linearly 
adding the additional spurious power as described above. 



4.5. Parametric approach to systematic effects 

In general, we do not use in this work a precise (still not completely available) description of the considered systematic 
effects, but only suitable representations of them. Therefore, we will use our estimations to understand if the considered 
classes of systematic effects may significantly affect the cosmological exploitation of Planck data with respect to the 
determination of cosmological parameters possibly by rescaling the estimation quoted above. This is done with the 
aim of identifying the corresponding levels at which it is necessary to control the systematic effects in order to avoid to 
spoil the scientific accuracy of the Planck data. We technically implement this rescaling with a multiplicative constant 
on the residuals of systematic effects on the C't described in the previous sections. 
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5. MOCK DATA PRODUCTION AND LIKELIHOODS 

We describe in this section the theoretical basis of our simulations when experimental errors are treated as statistical 
noise. This includes the instrumental white noise as well as the residuals from systematic errors and foregrounds as 
described in the previous section. In other words, we assume that the noise contribution to the observed CMB skies 
due to systematic errors can be precisely assessed, thanks to suitable procedures such as cleaning simulations in case 
of foreground residuals. 

Let us denote the fiducial theoretical multipoles as Cf T , Cj E , C EE , C BB and the (possibly ^—dependent) noise 
covariances as Nj J ', N EE and N BB . For instance, in the case of foreground residuals, we would have 

N xx ' = w e R xx ' + n xx ' , X, X' = T, E, B (9) 

where wi is the window function in multipole space, n xx is the white instrument noise and R xx are appropriate 
quantities which can be estimated through map cleaning simulations. In any case we assume here that N BB = N EE 
and that Nf E = Nf B = N EB = 0. 

Thus, the full covariances of the T — E fluctuations read 

WI CJ T +NJ T w e Cf E \ _ (C+ 0\ vt nm 

m CJ E Wl c EE + n ee ) =1Ze {o cj (10) 

where TZi are suitable rotation matrices and TZj stands for the trasposed matrix of TZp. 
The B fluctuations are decoupled, that is Cf E = Cf T = and have full covariance we Cf B + N BB . 

With these notations, a possible observed set of fluctuation amplitudes reads 

where g x m , X = +, — , B, are independent centered unit Gaussians, that is 

{oL) - , (gf m g x m ,) = 6 U > S mm , S xx ' (12) 
The amplitudes A x m are assumed real (which is always possible for integer weights £). 
The corresponding observed multipoles Ci (sometimes called pscudo-C^) arc 

m 

Consider now the multipoles C xx produced by some test cosmological model. In the approximation exploited in 

' 

this work (all-sky coverage and uniform sensitivity), the likelihood of such multipoles, given the observed C e , can 
be written as 

L = exp(-i X |, S -|xl) (14) 

where 



tr {C^Ct) - lo f 



detCe 
detCe 



(15) 



and Ci, Ci are the 2x2 matrices 



r [wi CJ T + NJ T Wl CJ E \ (cT C T i E \ n fi x 

^ \ Wl CJ E Wi C EE + N EE ) ' C * - [c™ cf E ) ( } 

while 

(jBB 

X * = ~ (~<BB I KfBB 

Wl u g + i\ e 

This setup is the minimal one but has the disadvantage that one cannot change the noise extraction while keeping 
the sky fixed. To allow such a possibility we consider independent signal and noise extractions as follows. First, we 
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introduce new rotation matrices TZi and eigenvalues Cf 1 so that now only the T — E signal covariances are diagonalizcd: 
Then, we double all Gaussian extractions by writing the fluctuations as 



Ala = V^~t «L + \JNJ T hj m , Af m = af m + y N EE hf m , Af m = ^ af m + y 1 N BB hf m , (19) 

where T, E and B stand for temperature. E polarization and B polarization, respectively, and the amplitudes o^ m are 
given by: 



a E 



tin 




B - . /„„„ nBB„B 



m Cf B gf m (20) 



The new independent centered unit Gaussians hX , Y = T, E, B, are independent from the previous set gf , X 
+,-,B. 



The pseudo-Cf , (that is Ci), can now be written as 



Cf = w e C xx + 2y/w e Q xx + N xx P xx , X = T,E,B 

C T t E = Wl Cj E + y/wt N EE Qj E + ^/w e QfT + VN TT N EE pTE 

where 



(21) 



with X, X' = T, E, B. 



This second approach allows to use the same sky and different noises, which is needed when cumulative channels 
are considered to reduce noise effects. For instance, a cumulative channel formed by the LFI at 70 GHz and the two 
HFI channels at 100 GHz and 143 GHz is obtained by simply summing the x 2 's of Eq. (|T5|) relative to these three 
channels. 

The above setup is based on the assumption that the noise contribution of systematic errors is precisely assessed. 
If this would not be the case, bias effects would be induced. This can be simulated in the likelihood x 2 's of Eq. (|15|) 

by using different noises in the C e and in the covariance built with the test multipoles C xx . That is, one should 

make (small) variations from Nf x to some N'f X in Eq. (fl6| while keeping them fixed in Eq. (Q~9|) (or viceversa) and 
study their impact on the parameter determination of the test cosmological model. 

6. FORECAST PRECISION OF THE PLANCK MEASUREMENTS FOR THE COSMOLOGICAL PARAMETERS WITHOUT 

SYSTEMATICS 

In our Monte Carlo Markov Chains (MCMC) simulations we take as fiducial model the ACDMr model, that is 
the standard ACDM model augmented by the tensor-to-scalar ratio r as described in the introduction. We consider 
MCMC simulations with both the ACDMr and the ACDMrT model. We denote by ACDMrT, the ACDMr model in 
which we impose the double-well inflaton potential given in Eq. (fTJ), as described in the introduction and in Sect. [2] 

We consider two sets of best fit fiducial values for our parameters, as listed in Tabled where also the values of few 
other derived parameters are shown for illustrative purposes. Since r = in the first set, the model is just the ACDM 
model. In the second set the values r = 0.0427 and n s = 0.9614 are chosen to lay on the theoretical curve r = r(n s ) 
dictated by the dou ble-well inflaton potential and they correspond to the best fit value y = 1.26 for the coupling 
(|Destri et al.ll2008al ) within the Ginsburg-Landau effective theory approach (see Eq. (p}). 

We then provide estimates of the errors in the measurements of the cosmological parameters in the following way: 

• We produce one sky (mock data) for Cf T , Cf E , Cf E and Cf B from the ACDM or the ACDMrT models (see 
Table [2]) according to the procedure described in Sect. [5j 

• We run Monte Carlo Markov Chains from this sky and obtain the marginalized likelihood distributions for the 
cosmological parameters (fib ft, 2 , il c h 2 , 9 r, 51a, Age of the Universe, z re , Hq, A s , n s and r) in the two 
models ACDMr and ACDMrT. To be precise, we study the independent ACDMr parameters with the mock 
data produced from ACDM (first row of Table [2]) and the independent parameters of both ACDMr and ACDMrT 
with the mock data produced from ACDMrT (second row in Table [2]). 
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We consider two choices for the Ci— likelihood, one without the B modes and one with the B modes and take into 
account the white noise sensitivity of Planck (LFI and HFI) in the 70, 100 and 143 GHz channels. We also consider 
a cumulative channel whose % 2 is the sum of the x 2 's of the three channels above. When using different channels in 
the MCMC analysis, we use different noise realizations while keeping the same sky, that is the same realization of the 
Gaussian process that generated the primordial fluctuations. 

In our MCMC analysis we always take standard flat priors for the cosmological parameters. In particular we assume 
the flat priors < r < 0.2 in the ACDMr model and < r < 8/60 according to the theoretical upper limit for r in 
the ACDMrT model. 

Our findings are summarized in Figs. [3]|6] where the marginalized likelihood distributions of the cosmological pa- 
rameters are plotted for several different setups. In Tables [2] to H] we list the corresponding relevant numerical values. 

Clearly, in the case of the ratio r, due to the specific form of its likelihood distribution it is more interesting to 
exhibit upper and lower bounds rather than mean value and standard deviation as in Tables [3] and 2] We report the 
upper bounds and, when present, the lower bounds in Tables [5] and [6] 

Notice that best and mean values reported here for r and the other cosmological parameters do not correspond to 
the real sky but to mock MCMC generated skies as explained above. However, the deviations between the best and 
the fiducial values are relevant indicators as well as the lower and upper bounds on r and the standard deviations. 

The fact that the fiducial and mean values of r are very close and that the standard deviation Ar max of the 
distribution of maximum values r max coincides with the mean value of the standard deviation of r indicate that Planck 
can provide measurements of high quality for r. 

As expected from the relative difference in sensitivity, typically the distributions obtained with the HFI-100 and 
HFI-143 channels agree very well while differing markedly from those obtained with the LFI-70 channel. Quite often 
the higher noise level in LFI-70 determines also shifts in the peak positions with respect to the other two channels. 
These shifts are within a la deviation in the LFI-70 distributions and represent therefore normal statistical fluctuations. 
The LFI-70 distribution on r when the fiducial value is r = 0.0427 does not exhibit a peak due to the sensitivity of 
this channel. 

As expected, whenever the probability distribution for a given parameter is close to Gaussian, the cumulative channel 
produces a distribution that is narrower than the narrowest distribution produced by any individual channel. This 
applies to all parameters in the ACDMr model, including r, which has a distribution close to a left -truncated Gaussian 
for both fiducial values used. 

In the ACDMrT model the relation between n s and r is non-linear and there are theoretical upper limits on n s and 
r (see Fig. [5]). These features introduce non-Gaussianities in the distributions and eventually also affect any other 
cosmological parameter having a sensitive correlation with n s and/or with r, such as fib h 2 , A s or some other derived 
parameters. Thus, the cumulative channel provides some parameter distributions which are larger than those of the 
HFI-143 channel, because it is affected by the LFI-70 channel, which is less sensible to constrain r well within its 
theoretical prior. This effect can be very well appreciated from Fig. [6] in which the likelihood distributions of the 
coupling constant y are plotted for the various cases considered. 

The very limited relevance of the B modes for the ACDMrT model is in principle expected because the n s value 
fixed by the T modes essentially determines r through the theoretical constraint. This property of the B modes shows 
up clearly from the figures and the tables. 

6.1. Forecasts of the Planck measurements with the ACDMr model 

The obtained best fits, mean values and standard deviations for the cosmological parameters are presented in Tables 
[3] and H] for the two cases considered: without B modes and with B modes included in the C;-likelihood respectively. 
Each table displays the values obtained for the two test models: ACDMr and ACDMrT. Two simulated skies (mock 
data) are considered: one with fiducial value r = (ACDM), and one with fiducial r = 0.0427 (ACDMrT). The 
complete sets of fiducial values used are given in Table [2j All values are rounded to order 10 -4 to the nearest value. 

From the Tables [3] to [6] we see that for the ACDMr model the best result, namely a peaked distribution for r, is 
obtained for a non-zero fiducial value for r (r = 0.0427) and with the B modes included in the likelihood. In this case, 
upper and lower bounds on r are obtained 

0.013 < r < 0.045 at 95%CL 

with the best values 

r = 0.0594 and n s = 0.9604 (without B modes); r = 0.0240 and n s = 0.9597 (with B modes). 

For a fiducial value r = 0, with or without B modes, the ACDMr distributions peak at the value r = 0, as seen from 
the upper panels of Figs. [5J31 Upper bounds on r are obtained in this case. They result to 95% CL 

r < 0.068 without B modes; r < 0.016 with B modes 

with the best values 
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r = 0.0041 and n s = 0.9549 (without B modes); r = 0.001 and n s = 0.9606 (with B modes). 

The results on n s practically do not change by including or not the B modes (compare the upper panels of Figs. [3] 
andE]). 

The upper bound on r and the best value of n s do not need the inclusion of B modes and can be obtained for fiducial 
r = 0. These values can be obtained and trusted without including the ACDMrT model, the pure ACDMr model is 
enough to obtain them. 

We see a substantial progress in the forecasted bounds for r wit h respect to the WMAP+LSS data set for which: 
r < 0.20 in the pure ACDMr model (jKomatsu et al.l [20091 12010D . For Planck, with the B modes included and a 
nonzero r-fiducial value, we get peaked distributions for r with a nonzero most probable value, lower bounds for r 
(and an improvement of the r upper bound). This is obtained by only using the ACDMr model alone without any 
input from the inflation model. We see now in the following subsection how these forecasts can be still considerably 
improved by using the ACDMrT model. 

6.2. Forecasts of the Planck measurements with the ACDMrT model 

With the ACDMrT model (namely when the double-well inflaton potential is imposed), with or without B modes 
included, well peaked distributions for r arc obtained together with upper and lower bounds and best r values. We 
get a considerable gain for r with respect to the pure ACDMr model, as can be seen from Tables [3] to [6] and Fig. [5j 

The fiducial value for r (r = 0.0427) is well reproduced by the peak of the ACDMrT distribution both with and 
without B modes. The ACDMrT distribution for r peaks at the non-vanishing value theoretically associated with the 
fiducial value of n s . We get (at 95% CL) 

0.030 < r < 0.113 (without B modes) and 0.030 < r < 0.114 (with B modes), 

with the best values and 95% CL errors 

r = 0.0463 ± 0.0231 and n s = 0.9625 ± 0.0035 (without B modes); 
r = 0.0405 ±0.0230 and n s = 0.9608 ± 0.0033 (with B modes). 

The results with the ACDMrT model practically do not change by including or not the B modes in the likelihood 
as can be seen from the figures (compare for instance the upper panels of Figs. [3] and 2]) and from the tables. This is 
so since the ACDMrT model intrinsically carries a non vanishing ratio prediction, which shows up in agreement with 
the obtained marginalized distributions even without the inclusion of B modes. 

6.3. Conclusion 

The upper bound on r and the best value of n s do not require to include the B modes in the likelihood, and can be 
obtained with the ACDMr model alone (see Fig. [3]). 

For a fiducial value r = 0, mock Planck data with the ACDMr model alone, with or without B modes in the 
likelihood, provide only upper bounds on r and most probable values for n 8 . 

The same conclusions are true for a non vanishing fiducial value (r = 0.0427) without B modes and the ACDMr 
model alone. 

The inclusion of B modes for a non vanishing fiducial value (r = 0.0427) allows peaked marginalized distributions 
for r with the ACDMr model alone and a lower bound for r (see Figs. 02 14] and [5|) . 

Lower bounds on r and most probable r values are always obtained (with or without the B modes) for the ACDMrT 
model (see Fig. [5]). 

The ACDMrT model provide in all the cases, with or without the B modes included, well peaked distributions for 
r on nonzero values r ~ 0.04. 

In summary, we find that the inclusion of the theoretical model greatly help the recovery of the r parameter. We 
also remark that the model is falsifiable in the case of constraints on n s and r not compatible with the banana shape 
of the considered framework. 

7. FORECAST PRECISION OF THE PLANCK MEASUREMENTS FOR THE COSMOLOGICAL PARAMETERS WITH TOY 

MODEL SYSTEMATICS 

7.1. Foreground residuals without bias 

We computed the cumulative marginalized likelihoods from the three channels including foreground residual^] for 
the cosmological parameters in the ACDMr and the ACDMrT models with B modes and fiducial ratios r = and 
r = 0.0427. 

12 Foreground residual exploited in this work overwhelms that coming from systematics earlier discussed. We included the latter in some 
representative tests, not reported here for sake of conciseness, finding that, as expected, it does not change significantly the conclusions 
derived taking into account the foreground residual only. 
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The foreground residuals arc introduced as an additional statistical error. We evaluate these statistical errors 
following the discussion and the toy model presented in Sect. [4] in which a worst case for considering the residuals is 
derived (as well as a best case, and a intermediate or middle case for the residuals). 

We plot the likelihoods for the cumulative of the three channels in four cases (see Figs. [7][8|): 

• (a) without residuals; 

• (b) with 30% of the toy model residuals in the TE and E modes displayed in Fig. [2] and 16 /J.K 2 in the T modes; 

• (c) with the toy model residuals in the TE and E modes displayed in Fig. [2] and 160 fiK 2 in the T modes; 

• (d) with 65% of the toy model residuals in the TE and E modes displayed in Fig. [2] and 88 \iK 2 in the T modes 
rugged by Gaussian fluctuations of 30% relative strength. 

The likelihood distributions with and without B modes result almost the same when including the residuals. Only 
the cosmological parameters sensitive to the B modes do appear to be affected by the residuals, namely, r, z re and r. 
This is so in the ACDMr model in which for a fiducial value r = 0, the upper bound in r does change by including 
the residuals (see Fig. [7]). This change is smaller for a fiducial ratio r = 0.0427. In this case, the presence of a lower 
bound for r (at 68% CL) does remain even by including the residuals. The lower bound remains at 95 % CL in the 
best case smooth residuals. The main numbers are displayed in Tables [JJ to [5] 

The r distribution for fiducial ratio r = 0.0427 is much clearly peaked on a value of r near the fiducial one in the 
ACDMr T model than in the ACDMr model (compare Figs. and 01. In any case, the best value for r in the presence 
of residuals is about r ~ 0.04 (near the fiducial value) both for the ACDMr and the ACDMrT models. 

The ACDMrT model turns to be robust, it is very stable with respect to the inclusion of the residuals as its 
marginalized likelihood distributions do not change (and we have seen in Sect. [5] that they do not change neither with 
respect to the inclusion or not of B modes). The main numbers are included in Tables [7] to GO We see in the ACDMrT 
model that we have again for r at 95 % CL: 

0.028 < r < 0.116 with the best value r = 0.04 . 

In summary, foreground residuals only affect B modes and therefore only the cosmological parameters sensitive to 
B modes are affected. 

7.2. Probability to detect r from B modes in the HFI-143 channel 

Figs. l9l and ITUl describe the probability of detection of r from B modes in the HFI-143 channel. 
In order to assess the probability for Planck to detect r we consider the B mode detection by the most sensitive 
HFI-143 channel. We follow the following method: 

• We extract 10 5 skies obtaining the corresponding multipoles Ai m from the ACDMrT model according to the 
procedure described in Sect. [5j We choose r = 0.0427 as fiducial value. 

• We compute all the corresponding likelihood profiles only for r. That is, freezing out all the other parameters 
to their fiducial values. 

• We compute the interesting properties of each likelihood profile, like the most likely value r max , the mean value 
fmeam the standard deviation Ar max of the r max distribution, the skewness and the kurtosis. This last measures 
the departure from a Gaussian likelihood. 

• We finally compute the 99% CL, 95% CL, and 68% CL lower bounds for r : rgg, rgs and r$s. 

In the five panels of Fig. [9] we plot the likelihood profiles for the different skies, r max and 6r max = r max — r mean . 
Notice that r max is an unbiased estimator of the true value, since its expectation value r mean throughout many skies 
coincides with the fiducial r. Ar max is the standard deviation of the r max distribution. We find that Ar max always 
agrees extremely well with the mean value of the standard deviation of r in each likelihood profile for each different 
sky. This fact means that asymptotically for a large number of skies the width of the r profile is an unbiased estimator 
of the actual uncertainty in r. 
All these results were obtained for a level of foreground residual equal to 30% of the toy model displayed in Fig. [2] 
We plot in Fig. ITOlthe 99% CL 95% CL, 68% CL lower bounds for r: rgg, r g5 and r 68 , respectively as functions 
of the fraction of foreground residual of the worst case. These lower bounds are consistent, since they fail more or 
less 99%, 95%, 68% of the sky extractions. This last property is true only if the prior r > is not enforced in the 
likelihood. That is why we get a non-zero likelihood on negative values of r. Of course, only positive values of r arc 
meaningful and this allows us to define the probability of detection of r, to 99% CL, 95% CL and 68% CL, as the 
fraction of skies which gives positive 99% CL, 95% CL and 68% CL lower bounds, respectively. 
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These probabilities of detection are displayed in Fig. [TUJ At the level of foreground residual equal to 30% of the 
considered toy model only a 68% CL (one sigma) detection is very likely. For a 95% CL detection (two sigmas) the 
level of foreground residual should be reduced to 10% or lower of the toy model displayed in Fig. [2] 

Lensing was not considered in the analysis on r— detection probability. Residuals and lensing affect the detection of 
B-modes in complementary ways, with the effect of residuals stronger than that of lensing. As a consequence, lensing 
plus residuals can spoil the detection of r even when residuals are assumed at the 30% level of the toy model displayed 
in Fig. [2] On the contrary, lensing in the absence of residuals still allows a detection of r. For example, several MCMC 
simulations show that our lower bounds on r are not significantly affected by lensing in the absence of foreground 
residuals. It should be clear that if the theoretical constraint r = r(n s ) of the ACDMrT model is imposed on the 
MCMC analysis, r has always well defined lower bounds regardless of lensing and/or residuals. 

The forecasted probability of detecting r is based on the statistics of the shape of the r -likelihood. This shape 
determines whether a detection of r can be claimed with a given confidence level. But real CMB experiments can 
observe only one realization (only one sample): the actual observed sky from which a single likelihood for the B-modc 
multipoles (and hence the value of r) is derived. So, the possibility of correctly inferring the value of r from one single 
(albeit very large) sample depends heavily on the sample itself, and therefore, in view of the detection probability 
found here, whether r will be or will be not detected depends also of a question on luck. 

In addition, the results for many skies presented in this section show the consistency of our whole approach to 
determine r. Similar results are valid for the other cosmological parameters. 

7.3. Bias effect in the foreground residuals 

Here we consider two extremal cases of bias, modelled as being an imprecise determination of the foreground residuals 
Re- We keep fixed the residuals introduced in the noise of the test covariances, while we change the residuals Rf x to 

some R'f X in the the noise of the observations, that is in Eq. Then, we write 

R'f x = R XX (1 + jflf) , X = T,E,B, (23) 
where for the numbers fie. we consider the two extremal cases: 

• (i) Independent flat random numbers fie from —0.5 to 0.5. Since in this case fie fluctuates randomly around zero 
the effect of the bias mostly cancel out and the cosmological parameters suffer little change as depicted in Fig. 

HU 

• (ii) Uniform ramps from — a x to a x as £ varies from 2 to 2100, a x varying randomly up to a 20% around 0.5 
with X = T, E, B. This means to choose 

^ = a (-L049-J- (24) 

Notice that there is always a non-zero value for a x despite its fluctuations and therefore there is a significant 
bias effect over the modes. In this case the bias depresses the estimated multipoles at low £ and increases them 
at high £. Thereby increasing the expected values of n s and depressing those of r. We see from Fig. [11] that r is 
not anymore detected in this case despite its fiducial value is always r = 0.0427. 

We consider a level of 30% of the toy model of foreground residuals displayed in Fig. [5] for the bias effect. In the case 
we change the overall sign of fi x in Eq. (|24[) . we introduce an additional spurious power in the estimation at low £ 
and a depression of the power at high £. This would erroneously increase the probability to detect r. 

The peaks in the cosmological parameters (with the exception of r) get shifted mainly due to the bias from the T 
modes. However, they stay within one or two a of the WMAP values. 

In case we do not add bias in the T modes, only r is affected significantly by the bias. Namely, in case (ii) without 
bias in T modes, all cosmological parameters except r peak practically at the same value as in absence of bias. On 
the contrary, the likelihood distribution for r is determined by the bias on the B modes and turns to be similar to the 
one in Fig. QT] for the bias case (ii). 

The bias introduced in case (ii) goes in the opposite direction to the theoretical double-well models where n s increases 
with r (see Fig. Q]). 

We only present here the bias for the ACDMr model. The likelihood distributions for the ACDMrT model including 
bias are similar to those of the ACDMr model except for r where a lower bound shows up due to the theoretical 
constraint. 

We only consider here two extreme cases of bias: case (i) where bias is practically harmless and case (ii) where it 
distorts significantly the cosmological parameters, especially r which is not anymore detected. 

8. FINAL CONCLUSION 

In this paper we provide a precise forecast for the Planck results on cosmological par ameters, in particular for th e 
tensor-to-scalar ratio r. These new forecasts go far beyond the published ones (see e.g. iPlanck Collaboration! (|2006ft . 
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iColombo et al.l (|2009D ) and pave the road for a promising scientific exploitation and interpretation of the Planck data 
(once cleaned from the different astrophysical foregrounds). 

We appropriately combined the following, as main ingredients: the current public available knowledge of Planck 
instrument sensitivity and a reason able toy model estimation of the residuals fro m systematic errors and foregrounds; 
the highly predictive theory setup ([Bovanovskv et aTl 120091 : iDestri c t al. 2008a. b) provided by the Ginzburg-Landau 
approach to inflation to produce and analyze the skies (mock data) which allows a decisive gain in the physical insight 
and data analysis; precise MCMC methods to produce the skies (mock data) and to analyze them. This turns into an 
improvement in the physical analysis, in particular for the ratio r. 

It must be also stressed that, in the considered framework, better measurements for n s will improve the predictions 
on r from the T, TE and E modes even if a secure detection of B modes will be still lacking. We remark also that 
the model is falsifiablc in the case of constraints on n s and r not compatible with the banana shape of the considered 
framework. 

The lower bounds and most probable value inferred from WMAP for r (r ~ 0.04) in the considered framework 
support the search for B mode polarization in Planck data and the future CMB B oriented polarization missions. 
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Frequency channel 


143 GHz 


100 GHz 


70 GHz 


Angular resolution [arcmin] 


7.1 


9.5 


13 


ST per FWHM 2 pixel [pK\ 


4.2 


4.8 


17.2 


<5Q, <5U per FWHM 2 pixel \pK] 


8.1 


7.7 


24.3 



TABLE 1 

Planck PERFORMANCE IN THE THREE FREQUENCY CHANNELS EXPLOITED IN THIS WORK. THE AVERAGE SENSITIVITY PER FWHM 2 
RESOLUTION ELEMENT (<5T, <5Q, <5U) IS GIVEN IN CMB TEMPERATURE UNITS (i.E. EQUIVALENT THERMODYNAMIC TEMPERATURE) FOR. 28 

MONTHS OF INTEGRATION, ALMOST CORRESPONDING TO FOUR SKY SURVEYS. 





n b h 2 


n c K 1 


6 


T 


log(10 iu A s ) 


ris 


r 


n A 


Ho 




ACDM 


0.0223 


0.1079 


1.0387 


0.0864 


3.0561 


0.9613 





0.7463 


71.628 


10.399 


ACDMrT 


0.0224 


0.1112 


1.0410 


0.0821 


3.0629 


0.9615 


0.0427 


0.7364 


71.228 


10.062 



TABLE 2 

Fiducial parameters for the two considered models: the standard ACDM model and the ACDMrT model. In the ACDMrT 

MODEL WE CONSTRAIN r = r(n s ) BY THE DOUBLE-WELL INFLATON POTENTIAL GIVEN IN EQ. {TJ AS DEPICTED ON THE UPPER BORDER OF 

THE BANANA-SHAPED REGION FlG. Q] 



Sky data 


ACDM: fiducial r = 


ACDMrT: fiducial r = 0.0427 


test model 


ACDMr 


ACDMr 


ACDMrT 


without B 
modes 


best 


mean 


stddev 


best 


mean 


stddev 


best 


mean 


stddev 


fl b h 2 


0.0223 


0.0223 


0.0001 


0.0225 


0.0225 


0.0001 


0.0227 


0.0226 


0.0001 


n c h 2 


0.1085 


0.1082 


0.0007 


0.1118 


0.1116 


0.0007 


0.1112 


0.1114 


0.0010 


e 


1.0389 


1.0388 


0.0002 


1.0411 


1.0412 


0.0002 


1.0415 


1.0413 


0.0003 


T 


0.0832 


0.0833 


0.0027 


0.0858 


0.0857 


0.0027 


0.0871 


0.0872 


0.0028 


log[10 iu A s ] 


3.0479 


3.0618 


0.0054 


3.0703 


3.0697 


0.0056 


3.0722 


3.0729 


0.0055 


n s 


0.9549 


0.9609 


0.0021 


0.9604 


0.9607 


0.0022 


0.9625 


0.9619 


0.0035 


r 


0.0041 


0.0284 


0.0206 


0.0594 


0.0532 


0.0277 


0.0463 


0.0510 


0.0231 


n A h 2 


0.7436 


0.7451 


0.0038 


0.7344 


0.7351 


0.0040 


0.7388 


0.7369 


0.0059 




71.437 


71.577 


0.3610 


71.111 


71.165 


0.3638 


71.591 


71.397 


0.5288 


2tc 


10.128 


10.121 


0.2267 


10.375 


10.358 


0.2328 


10.414 


10.446 


0.2276 



TABLE 3 

Best FITS, MEAN VALUES and STANDARD DEVIATIONS for cosmological PARAMETERS WHEN B MODES ARE NOT INCLUDED IN THE 
C^-LIKELIHOOD. ALL VALUES ARE ROUNDED TO ORDER 10 -4 TO THE NEAREST VALUE AND CORRESPOND TO THE CUMULATIVE CHANNEL 
WHOSE X 2 IS THE SUM OF THE X 2 'S OF THE THREE CHANNELS HFI-100, HFI-143 AND LFI-70. 



Sky data 


ACDM: fiducial r = 


ACDMrT: fiducial r = 0.0427 


test model 


ACDMr 


ACDMr 


ACDMrT 


with B modes 


best 


mean 


stddev 


best 


mean 


stddev 


best 


mean 


stddev 


fl b h 2 


0.0223 


0.0223 


0.0001 


0.0225 


0.0225 


0.0001 


0.0226 


0.0226 


0.0001 


fl c h 2 


0.1081 


0.1084 


0.0007 


0.1118 


0.1118 


0.0007 


0.1117 


0.1114 


0.0010 


e 


1.0389 


1.0389 


0.0002 


1.0412 


1.0411 


0.0002 


1.0412 


1.0413 


0.0003 


T 


0.0834 


0.0840 


0.0027 


0.0866 


0.0857 


0.0027 


0.0865 


0.0873 


0.0028 


log[10 10 A s ] 


3.0620 


3.0624 


0.0053 


3.0721 


3.0703 


0.0054 


3.0720 


3.0728 


0.0056 


n s 


0.9606 


0.9603 


0.0021 


0.9597 


0.9602 


0.0021 


0.9608 


0.9621 


0.0033 


r 


0.0010 


0.0060 


0.0050 


0.0240 


0.0275 


0.0096 


0.0405 


0.0516 


0.0230 


fl A h 2 


0.7456 


0.7445 


0.0038 


0.7342 


0.7340 


0.0040 


0.7356 


0.7373 


0.0052 


Ho 


71.624 


71.521 


0.3556 


71.089 


71.075 


0.3570 


71.265 


71.430 


0.4931 


2tc 


10.126 


10.132 


0.2250 


10.443 


10.366 


0.2285 


10.397 


10.448 


0.2319 



TABLE 4 

Best FITS, MEAN values AND STANDARD DEVIATIONS FOR COSMOLOGICAL parameters WHEN B MODES are INCLUDED IN THE 
C^-LIKELIHOOD. ALL VALUES ARE ROUNDED TO ORDER 10~ 4 TO THE NEAREST VALUE AND CORRESPOND TO THE CUMULATIVE CHANNEL 
WHOSE X 2 IS THE SUM OF THE X 2 'S OF THE THREE CHANNELS HFI-100, HFI-143 AND LFI-70. 
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sky data 


ACDM: fiducial r = 


ACDMrT: fiducial r = 0.0427 


test model 


ACDMr 


ACDMr 


ACDMrT 




68% CL 


95% CL 


68% CL 


95% CL 


68% CL 


95% CL 


without 
B modes 


LFI-70 


r < 0.2 


r < 0.2 


r < 0.2 


r < 0.2 


r < 0.102 


r < 0.128 


HFI-100 


r < 0.068 


r < 0.124 


r < 0.097 


r < 0.155 


r < 0.047 


r < 0.111 


HFI-143 


r < 0.070 


r < 0.117 


r < 0.108 


r < 0.158 


r < 0.042 


r < 0.061 


Cumulative 


r < 0.036 


r < 0.068 


r < 0.066 


r < 0.102 


r < 0.051 


r < 0.113 


with 
B modes 


LFI-70 


r < 0.074 


r < 0.151 


r < 0.075 


r < 0.144 


r < 0.117 


r < 0.131 


HFI-100 


r < 0.012 


r < 0.029 


r < 0.037 


r < 0.065 


r < 0.049 


r < 0.112 


HFI-143 


r < 0.008 


r < 0.020 


r < 0.041 


r < 0.064 


r < 0.042 


r < 0.062 


Cumulative 


r < 0.008 


r < 0.016 


r < 0.032 


r < 0.045 


r < 0.052 


r < 0.114 



TABLE 5 

Upper bounds on r with all figures rounded upward to order 10~ 3 . Notice that the bound r < 0.2 is just the assumed 

PRIOR, WHICH GETS SATURATED BY THE ACDMr TEST MODEL IN THE LFI-70 CHANNEL WHEN B MODES ARE ABSENT. THE LIMITS IN THE 
CASE OF THE ACDMrT TEST MODEL WITH FIDUCIAL r = 0.0427 IN THE LFI-70 CHANNEL ARE NOT REALLY SIGNIFICANT IN VIEW OF THE 

SHAPE OF THE CORRESPONDING LIKELIHOOD DISTRIBUTION (SEE FlG. . 





sky data 


ACDMrT: fiducial r = 0.0427 


test model 


ACDMr 


ACDMrT 




68% CL 


95% CL 


68% CL 


95% CL 


without 
B modes 


LFI-70 






r > 0.067 


r > 0.035 


HFI-100 


r > 0.051 




r > 0.034 


r > 0.024 


HFI-143 


r > 0.034 


r > 0.026 


r > 0.030 


r > 0.025 


Cumulative 


r > 0.024 




r > 0.034 


r > 0.030 


with 
B modes 


LFI-70 






r > 0.046 


r > 0.024 


HFI-100 


r > 0.020 




r > 0.034 


r > 0.025 


HFI-143 


r > 0.026 


r > 0.013 


r > 0.033 


r > 0.026 


Cumulative 


r > 0.022 


r > 0.013 


r > 0.039 


r > 0.030 



TABLE 6 

Lower bounds on r with all figures rounded downward to order 10~ 3 . These bounds are assumed significant only when 

THE LIKELIHOOD AT r = IS LESS THAN exp( — 1/2) = 0.6065 ... OF ITS MAXIMUM FOR 68% CL OR LESS THAN exp( — 1) = 0.3678 . . . FOR 

95% CL. The limits in the case of the ACDMrT test model with fiducial r = 0.0427 in the LFI-70 channel are not really 

SIGNIFICANT IN VIEW OF THE SHAPE OF THE CORRESPONDING LIKELIHOOD DISTRIBUTION AS CAN BE SEEN FROM FlG. [5] In THE ACDMr 

model, the entries left empty in the table correspond to the cases where there are no lower bounds on r (as can be seen 

from Figs. [3]and[4|. 



sky data 


ACDM: fiducial r = 


ACDMrT: fiducial r = 0.0427 


test model 


ACDMr 


ACDMr 


ACDMrT 




best 


mean 


stddev 


best 


mean 


stddev 


best 


mean 


stddev 


no residuals 


0.0010 


0.0060 


0.0050 


0.0240 


0.0275 


0.0096 


0.0405 


0.0516 


0.0230 


best case smooth 


0.0040 


0.0222 


0.0160 


0.0448 


0.0504 


0.0238 


0.0465 


0.0516 


0.0235 


middle case rugged 


0.0024 


0.0230 


0.0188 


0.0431 


0.0472 


0.0261 


0.0344 


0.0513 


0.0234 


worst case smooth 


0.0083 


0.0250 


0.0160 


0.0436 


0.0480 


0.0275 


0.0387 


0.0518 


0.0250 



TABLE 7 

Best fits, mean values and standard DEVIATIONS of the ratio r WHEN B MODES are included in THE CUMULATIVE — LIKELIHOOD 
AND FOREGROUND RESIDUALS ARE TAKEN INTO ACCOUNT. ALL FIGURES ARE ROUNDED TO ORDER 10 -4 TO THE NEAREST VALUE. 



sky data 


ACDM: fiducial r = 


ACDMrT: fiducial r = 0.0427 


test model 


ACDMr 


ACDMr 


ACDMrT 




68% CL 


95% CL 


68% CL 


95% CL 


68% CL 


95% CL 


no residuals 


r < 0.008 


r < 0.016 


r < 0.032 


r < 0.045 


r < 0.052 


r < 0.114 


best case smooth 


r < 0.028 


r < 0.053 


r < 0.062 


r < 0.091 


r < 0.052 


r < 0.115 


middle case rugged 


r < 0.029 


r < 0.058 


r < 0.059 


r < 0.094 


r < 0.052 


r < 0.115 


worst case smooth 


r < 0.032 


r < 0.062 


r < 0.060 


r < 0.097 


r < 0.052 


r < 0.116 



TABLE 8 

Upper bounds on r when foreground residuals are considered, using the cumulative likelihoods and including B modes. 

All figures are rounded upward to order 10 -3 . 
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sky data 


ACDM: fiducial r = 


ACDMrT: fiducial r = 0.0427 


test model 


ACDMr 


ACDMr 


ACDMrT 




68% CL 


95% CL 


68% CL 


95% CL 


68% CL 


95% CL 


no residuals 






r > 0.022 


r > 0.013 


r > 0.039 


r > 0.030 


best case smooth 






r > 0.037 


r > 0.012 


r > 0.038 


r > 0.029 


middle case rugged 






r > 0.032 




r > 0.038 


r > 0.029 


worst case smooth 






r > 0.032 




r > 0.037 


r > 0.028 



TABLE 9 

Lower bounds on r when foreground residuals are considered, using the cumulative likelihoods and including B modes. 
These bounds are assumed significant only when the likelihood at r = is less than exp(— 1/2) = 0.6065 ... of its maximum 
for 68% CL or less than exp(— 1) = 0.3678. . . for 95% CL. Otherwise the entry is left empty in this table. All figures are 

rounded downward to order 10 -3 . 
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0.08 




0.06 



Fig. 1. — The universal banana region B in the (n s ,r)-plane setting N = 60. The upper border of the region B corresponds to the fourth 
order double- well potential expressed by Eq. JTJ. The lower border is described by the potential V(ip) = i«i 2 <p 2 *) for </? 2 < m 2 /A 



and V((p) = co for i/? 2 > m 2 /A l|Destri et alj|2009h . We display in the vertical full line the ACDMr value n s 
WMAP+BAO+SN data set. The broken vertical lines delimit the ±lcr region. 



= 0.968 ± 0.015 using the 
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Fig. 2. — Comparison between our model of foreground residual in the E mode (solid line) and in the TE mode (dashes) with a typical 
CMB angular power spectrum for the same modes and our fiducial B mode (i.e. with r = 0.0427) including (upper dashed-dot line) or not 
(lower dashed-dot line) the contribution by lensing. Foreground residuals in the B modes are assumed equal to those in the E modes. 
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ACDMr without B modes and fiducial r - 




log(10 10 A s ) 



ACDMr without B modes and with fiducial r - 0.0427 




log(10 10 A s ) 



Fig. 3. — Marginalized likelihood distributions, without including B modes, of the cosmological parameters for the ACDMr model. We 
display the distributions for each of the three channels HFI-100, HFI-143 and LFI-70 and for the cumulative of the three channels. The 
fiducial values are indicated by a vertical thin black line. The fiducial value for the ratio is r = in the upper panel and r = 0.0427 in the 
lower panel. Notice that the latter fiducial value is smaller than the peaks of the marginalized distribution. This is due just to statistical 
fluctuations since we are considering only one sky. The upper panel figures imply the upper bound r < 0.068 at 95% CL without B modes. 
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Fig. 4. — Marginalized likelihood distributions, including B modes, of the cosmological parameters for the ACDMr model. We display 
the distributions for each of the three channels and for the cumulative of the three channels. The fiducial values are indicated by a vertical 
thin black line. The fiducial value for the ratio is r = in the upper panel and r = 0.0427 in the lower panel. Notice that the latter fiducial 
value is larger than the peaks of the marginalized distribution. This is due just to statistical fluctuations since we are considering only one 
sky. The lower panel figures give upper as well as lower bounds for r : 0.013 < r < 0.045 at 95% CL. The results on n s practically do not 
change by including or not the B modes as we see from the upper panels of Figs. [3] and [4] 
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ACDMrT without B-modes and with fiducial r = 0.0427 
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ACDMrT with B-modes and fiducial r = 0.0427 




LFI-70 

HFI-100 

HFI-143 

Cumulative 

3.02 3.06 3.1 3.14 0.95 0.96 0.970 0.05 0.1 0.15 
log(10 10 A s ) "s ' 

Fig. 5. — Marginalized likelihood distributions of the cosmological parameters for the ACDMrT model in which the double— well inflation 
theoretical model is imposed. The MCMC analysis includes (does not include) B modes in the lower (upper) panel. We display the 
distributions for each of the three channels and for the cumulative of the three channels. The fiducial values are indicated by a vertical 
thin black line. The fiducial value of r is well reproduced by the peak of the ACDMrT distribution in the case of the HFI-100, HFI-143 
and cumulative of the three channels both with and without B modes. Considerable gain is obtained with respect to the ACDMr model. 
Upper and lower panels show quite similar results showing the stability of the ACDMrT model with respect to the inclusion of the B 
modes. Comparison with the lower panel of Fig. [4] shows that considerable gain is obtained with respect to the ACDMr model. 
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without B-modes 



with B-modes 




Fig. 6. — Marginalized likelihood distributions of the coupling constant y of the double— well quartic inflaton potential in the ACDMrT 
model. The MCMC analysis includes (does not include) B modes in the right (left) panel. We display the distributions for each of the 
three channels and for the cumulative of the three channels. The fiducial values are indicated by a vertical thin black line. The fiducial 
value of y is relatively well reproduced by the peak of the distribution in the case of the HFI-100, HFI-143 and cumulative channels. 
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ACDMr with B-modes, fiducial r = and foreground residuals 
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Fig. 7. — Cumulative 3— channel marginalized likelihood distributions, including B modes and foreground residuals, of the cosmological 
parameters for the ACDMr model. The fiducial ratio is r = in the upper panel and r = 0.0427 in the lower. We plot the distributions in 
four cases: (a) without residuals, (b) with 30% of the toy model residuals in the TE and E modes displayed in Fig. [2] and 16 fiK 2 in the T 
modes, (c) with the toy model residuals in the TE and E modes displayed in Fig. [2] and IGOfiK 2 in the T modes, (d) with 65% of the toy 
model residuals in the TE and E modes displayed in Fig. [2] and 88(iK 2 in the T modes rugged by Gaussian fluctuations of 30% relative 
strength. 
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Fig. 8. — Cumulative marginalized likelihoods from the three channels for the cosmological parameters for the ACDMrT model including 
B modes and fiducial ratio r = 0.0427 and the foreground residuals. We plot the cumulative likelihoods in four cases: (a) without residuals, 
(b) with 0.3 of the worst case residuals in the TE and E modes and 16fiK 2 in the T modes, (c) with the worst case residuals in the TE 
and E modes and 160/iif 2 in the T modes, (d) with 65% of the toy model residuals in the TE and E modes displayed in Fig. [2]and 88fj,K 2 
in the T modes rugged by Gaussian fluctuations of 30% relative strength. 
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Fig. 9. — Upper left panel: the likelihood profiles for the different skies. Upper right panel: r max and Sr max = r max — r mean and the 
fiducial r, r f id- Lower left panel: the skewness and the kurtosis. Lower middle panel: the 99% CL, 95% CL, and 68% CL lower bounds for 
r : rgg, rgs and r68- Lower right panel: the standard deviation std of the r distributions for each sky and the standard deviation Ar m ax 
of the r max distribution. These results correspond to a level of foreground residual equal to 30% of the toy model displayed in Fig. [2] The 
sensitivity of the 143 GHz channel is exploited here. 
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Fig. 10. — 99% CL, 95% CL and 68% CL lower bounds for r as functions of the fraction of foreground residual of the worst case. For 
30% foreground residual case only a 68% CL detection is very likely. For a 95% CL detection the level of foreground residual should be 
reduced to 10% or lower of the toy model displayed in Fig. [2] The sensitivity of the 143 GHz channel is exploited here. 



ACDMr with B-modes, fiducial r = 0.0427, bias and foreground residuals 











// 


V * 










II 




1 


/V. 


Jr 










it 
ii 


' V 


/' ' 










0.02240.02260.0228 

a tr 

b 



0.11 0.1,12 0.1141.0406 1.041 1.0414 1.0418 0.07 0.08 0.09 0.1 

£2, tr e x 








' * / 




i \ / 


\\ 




\\ 


/ / \ 


x\ 


'AT s 


\ 



0.72 0.73 0.74 0.75 13.65 13.7 13.75 













\ 


III 




7' 
J' 


* V 


9 10 11 



a. 



Age/GYr 




3.04 3.06 lr 3.08 3.1 0.955 0.96 0.965 0.97 
log(10 10 AJ n 





best case smooth 
b.c.s. with bias (i) 
b.c.s. with bias (ii) 
worst case smooth 



Fig. 11. — Likelihoods profiles with bias and foreground residuals including B modes for the ACDMr model and fiducial value r = 0.0427. 
We plot the best and worst smooth case s of the residuals without bias and the best smooth case for the residuals including the bias 
according to the cases (i) and (ii) in Sect. 17.31 namely small and large bias cases respectively. 



